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Abstract 


We study a two-dimensional semi-totalistic binary cell-state cellular au- 
tomaton, which imitates a reversible precipitation in an abstract chemical 
medium. The systems exhibits a non-trivial growth and nucleation. We 
demonstrate how basic computational operation can be realized in the 
system when the propagation of the growing patterns is self-restricted by 
stationary localizations. We show that precipitating patterns of different 
morphology compete between each other and thus implement serial and 
non-serial logical gates. 


Keywords: cellular automata, chemical reaction, growth, nucleation, waves, 
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1 Introduction 


Constraining the media geometrically is a common technique used when de- 
signing computational schemes (particularly implementing logical circuits) in 
spatially extended non-linear media. For example ‘Strips’ or ‘channels’ are 
constructed within the medium (e.g. excitable medium) and connected to- 
gether, typically using arrangements such as T-junctions. Fronts of propagating 
phase (excitation) or diffusive waves represent signals, or values of logical vari- 
ables. When fronts interact at the junctions some fronts annihilate or new 
fronts emerge. They represent a result of the computation. The approach is far 
from being elegant (compared for instance to the collision-based paradigm [[]) 
but simple and practical enough to be used by experimental scientists and en- 
gineers. This method of constraining the excitable media geometrically has 
enabled the design of a range of computing devices in the chemical laboratory. 


This has included logical gates [23], diodes [17] [9] 20], counters [TI], coin- 
cidence detectors [I2] and memory [19] all carried out using analogues of the 
Belousov-Zhabotinsky (BZ) medium. 

Until very recently precipitating chemical systems had not been used in 
unconventional computing. However, previously we demonstrated in chemical 
laboratory experiments [3] that an XOR gate can be constructed in the reaction 
between palladium chloride and Potassium iodide (the palladium processor) via 
the use of a simple T-junction constructed using a substrate loaded gel. In the 
reaction precipitation fronts stop propagating when they encounter each other 
leaving a precipitate free zone between the competiting fronts (competition for 
the substrate). The ideas was further developed in the construction of logical 
gates in propagating slime mold [26] (repulsion of growing patterns). 

As discussed in [4] there are certain (dis)advantages of excitable and pre- 
cipitating chemical media. Excitable systems are always (assuming unlimited 
supply of resources and removal of by-products) returning to their original rest- 
ing state, and thus are capable of processing multiple streams of information. 
However they do not have static memory. Precipitating systems are disposable 
(it takes too much energy to reverse precipitation). However precipitating sys- 
tems always maintain the result of the computation in their long-term memory. 
So, both types of systems are equally ‘useful’. 

In the present paper we investigate the possibility of computation in self- 
constrained reaction-diffusion systems. We develop and analyze a cellular au- 
tomaton model where the channels are constructed using stationary localizations 
and computation is implemented by propagating patterns of precipitation. 


2 The cellular automaton 


Previously we have proposed a detailed classification of two-dimensional semi- 
totalistic cellular automaton functions [5], which imitates the processes of dis- 
sociation and association in spatially extended nonlinear media. Every cell of 
the automaton array has eight neighbours and takes two states, 0 and 1. Cells 
update their states simultaneously in discrete time, depending on the sum of 
the states of their neighbours, and by the same cell-state transition rule. A 
cell in state 0 takes state 1 if the number of neighbours in state 1 belongs to 
the interval [01,02], otherwise the cell remains in the state 0. A cell in state 1 
continue in state 1 if the number of its neighbours in state 1 belongs to interval 
(51, 62], otherwise the cell changes the state 0. The transition parameters satisfy 
the following condition: 0 < 6; < 62 < 8 and 0 < 6, < 69 < 8. This way, the 
rules can be written as 0106169. 

Such a cell-state transition function can be interpreted as a simple discrete 
model of a quasi-chemical system with substrate ‘0’ and reagent ‘1’, where 
[01,02] is an interval of reaction (association or precipitation), and [ô4, 9] is an 
interval of dissociation. These family of rules includes Conway’s Game of Life 
(GoL), when intervals [6), 52] and [61,42] are interpreted as intervals of survival 


Sabena 


Figure 1: Patterns constructed with the evolution rule 2422. First case show an 
initial random perturbation into a box of 10 x 10 cells, initial density of %50, 
evolution in 335 times with 105,456 live cells. Second case show reactions of two 
gliders with a small oscillator period two, evolution in 360 times with 107,519 
live cells. 


and birth respectively[] 

We have demonstrated that this family of cellular automata is capable of 
constructing Voronoi diagrams [4], information transfer in the form of mobile 
self-localizations [15] and self-development of sophisticated patterns from mini- 
mal seed conditions [5]. 

In this way, we have also discovered and verified in computational exper- 
iments that a bigger cluster of rules Life dc22, when d and c takes values of 
2 to 8 and d < c uniquely supports mobile self-localizations (gliders) ?| almost 
frequently growing in square patterns, and complex patterns combining highly- 
disordered and ordered zones they are reached by initial perturbations or from 
reactions of mobile self-localizations [5]; see some examples in Fig. 

An special attention was concentrated to develop computations in the do- 
main of GoL with the also called “Diffusion Rule” [5]. However common self- 
localizations were also involved into the big cluster dc22. In this way, a subset 
of rules called Life 2c22 was analyzed independently founding mobile and sta- 
tionary self-localizations although they were very limited to project possible 
computations like in GoL. So, a significant problem was stop supernova ex- 
plosions, but like was reported in [I6] an special pattern well-known as “still 
life” pattern in GoL, exist in 2c22. Generally, any CA that grown quickly was 
very complicated to stop. In this paper, we report as stop these explosions 
constructing “walls” in the rules 2c22 when c > 3. 


Figure 2: Stationary self-localizations containing unlimited growth, two cases 
they are displayed. First case show an static universe that conserve its form 
even while a chaotic growth is present. Second case show an empty universe 
where a chaotic growth is evolving out of this area. 


Thus a stationary self-localization is defined into a square of 6 x 6 cells. A 
wall is defined like the volume of this stationary structure plus a cell, i.e., volume 


1Life Info http: //www.pentadecathlon.com/ 


‘http: //uncomp.uwe.ac.uk/genaro/diffusionLife/life_2c22. html 


+1 cell. This construction permit contain any information or perturbation as 
an static universe (see Fig. |2). 

Stability with unexcited states in initial conditions with excited states, in 
general was related to nucleation phenomenon [I3], probabilistic analysis with 
mean field theory [14] [I8] confirm that these rules 2c22 when c > 3 most concen- 
trations of cells in state 1 domain, gradually its curve growth fixing stable fixed 
points in state 1. But, very small initial densities in state 1 (produced by un- 
stable fixed points) produce a nucleation phenomenon from small perturbations 
or with reactions of self-localizations. Of course, it is related to minimum noise 
necessary into the evolution space to produce nucleation like in an excitable 
chemical reaction [8]. 


Figure 3: Enlarged photo of spiral wave-fronts in precipitating sytems: copper 
chloride and potassium hydroxide. Size is 1 x 1mm Øl. 


Complex growing patterns generated by rule 2c22 do reflect the surprising 
complexity of some natural precipitating systems, where the apparent spatio- 
temporal dynamics is morphologically comparable to that of excitable chemical 
systems. A remarkable example is shown in Fig. |3} There we can see complex 
structures formed when certain concentrations of a simple salt, copper chloride, 
immobilized in thin gel sheets is reacted with a solution of potassium chloride [f]. 

In the present paper we construct basic computing primitives using configu- 


rations of cellular automata governed by rules 2c22. We use the propagation of 
waves in channels of information with stationary and mobile self-localizations. 
Also, we will see how processing serial and non-serial logic gates. In simulations 
we have employed Golly system|}| 


3 Computing in rules 2c22 


To implement a universal computation, in terms of Boolean logical gates, in a 
spatially extended medium one must quantize information in the local propagat- 
ing patterns. To enable quantization we must prevent any initial configuration 
from growing unlimitedly. 

In automata rules 2c22 we have selected the rule 2622 to made our con- 
structions (clarity in their patterns). Of course, we could use their other four 
rules but it is not necessary because the other rules were just an adjustment 
to manipulate the same patterns, potentially all they are candidates to process 
the same operations and therefore proof its respective logic universality in 2c22 
when c > 3. 


(a) (b) 


Figure 4: Two types of patterns propagation in the information channels (bot- 
tom) and their initial configurations stimulated by a mobile self-localization 


(top). 


Thus we can make an impenetrable barrier of stationary self-localizations, 
so-called still life configurations as we saw in previous one section. Primitively, 
each stationary localization is an intersection of four segments of 1 state, each 
segment occupies six cells in state 1; the segments intersect so to enclose the 
2 x 2 domain of cells in state 0. Eventually, a family of patterns there is manip- 
ulating its width of channel and positions of mobile self-localizations. But, we 
have used minimum distance by mobile and stationary self-localization and it 
was calculated as: interval = volume still life — volume glider. Minimal space 
between walls limit number of patterns into of them. See exact configurations 
of the barrier localizations in Fig. 


http://golly.sourceforge.net/ 


Signals are represented by two types of patterns propagating along the chan- 
nels (Fig. |4| bottom). Both patterns are initiated by the same localization 
(Fig.|4} top), consisting of four cells in state 1, but differently positioned. Logi- 
cal value 0 is represented by regular wave fronts, which are generated by placing 
localization at the bisector of the channel’s longitudinal walls (Fig. (4p). To gen- 
erate less regular patterns (representing value 1) we shift the localization one 
cell away from the bisector (Fig. [4p). 

Truth values of Boolean variables are represented by these two propagating 
patterns. Computation occurs when patterns propagating in channels interact 
at the junctions of the channels. Configurations of cellular automata implement- 
ing the logical gates AND and OR are shown in Fig.|5| Automaton configurations 
before computation are shown on the left (initial state), besides you can see lo- 
calizations, representing logical values 0 and 1 depending its position or phase. 
The configuration at the end of computation, after the patterns have completed 
their propagation and interaction are shown on the right in both cases. 


HHRHHRHHHHHH 


gnnnnnnnnnHH 
Be s 
snnnnnNREE 


Re 
HHRHHHHHHH H 
wR 


# 
bi ghidetatas 4 
HHRHHHHHH 
# 


# 
LEsaReeee # 


© # 
HHRHHHHHHHHHH 


1/0 


# 
# 
# 
# 


RRR 
# 
# 
# 
# 
# 
# 
# 
# 
aR 


# 
gnnnnnnnnt 
nnnnnRRRERRH 


1/0 


HHRHHHHHHHHHH 


ž nnnnnRNRNRHH 
nnnnnRRBEN H 
Re 


nnnnnRBBEN H 
RE 


1/0 


# 
HHRHHHHHH 
HHRBHHHHH 
# 


# 
HHHHHHHH 

# 
nnnnnRRH 
# 


nunHH 


# 
# 
# 
# 
# 


ae 
HHRHHHHHHHH H 


ae 
HHHHHHHEHH H 
Sista teases st tet t 


Peereree eer t 


11% 1⁄1 


gnnnnnnnnnNE 1 nnnnnnnnRRNH 1 
HrnnERRENH H nnnnnnHHHH H 

HOR HOR 

HOR HOR 

g Bnaenene H Rete 

# 

H penne B nnnnRRHE 

HOt HR 

HOR HOR 
BRRHERNHEH H nnnnnRHBHE H 
© # Be # 1/0 
nnnnnnBRERRH nnnnnRRBERHH 
gnnnnnnnnnNE nnnnnnnBRRRH 1 
Btnnaseenn H nnnnnRRBEN H 

HOR OR 

HOR HOR 

# sunnnnng g sunnnnng 

# 

H anaes H panna 

# OF HOR 

OR HOR 
BERBHRNHHE H BERBERNHHE H 
He # Be # 1⁄1 
nennnRRRERHH nennnRRRRRRH 

(a) (b) 


Figure 5: Snapshots of cellular automaton configurations, corresponding to im- 
plementation of logical gates (a) AND and (b) OR. 


The channels and junctions, made of stationary localizations, are exactly 
the same for AND and OR gates? The difference is in the exact position of the 
localization-stimulus representing value 1. Compare configurations correspond- 
ing to input (01) in Fig. [6pb. To implement AND gate the localization repre- 
senting 1 must be shifted toward the localization representing 0 (shifted North) 
in Fig. [6h. While for the OR gate the localization representing 1 must shifted 
away from the localization representing 0 (shifted South). Such shifts lead to 
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Figure 6: Cellular automaton implementation of MAJORITY gate: (a) all cases 
corresponding to majority input value I/0, (b) all cases corresponding to ma- 
jority input value I/1. 


different propagating patterns, representing 1, which either lose (AND gate) or 
win (OR gate) in competition with the propagating pattern, which represents 0. 

Implementation of MAJORITY gate is shown in Fig. [6] The gate has three 
inputs: North, West and South channels, and one output: East channel. This 
gate was constructed when their three propagating patterns hit in the center of 
the channels thus the final result (O/) is the result of this collisions of waves. 
Some similarity to construct this gates can be founded with quantum-dot cellular 
automata [21]. 


4 Discussion 


We presented a cellular automaton model of an abstract quasi-chemical system, 
where selective and reversible precipitation can occur depending on the local 
concentration of the species. When the system in a quiescent state is locally 
disturbed, the disturbance grows unlimitedly. However, if it is surrounded by 
stationary localizations the disturbance propagates in a controllable manner. 
We demonstrate that such a property of a self-imposed geometrical constraint 
can be used to implement logical gates, and consequently circuits. The results 
of the paper give yet more evidence that precipitating chemical systems can be 
used to construct universal logical circuits. 

Self-restriction of the — otherwise — unlimited growth is yet another novel 
feature of the paper. Previous results, e.g. Banks cellular automata [6] or sand 
pile automata [10] have already exploited stationary perturbations, or wires, to 
guide traveling patterns/signals. However, in their work, signals were rather 
similar to electrical potential propagating along a wire of conductive material. 
This way only the rule 2622 demonstrates computation in the system with un- 
limited growth self-constrained by stationary localizations. However, the cluster 
of rules 2c22 (c > 3) have the same capacities to support logic universality. 

The computation in the 2c22 medium is based on competition between two 
types of propagating patterns. We believe the theoretical results discussed will 
demonstrate their practical value in experimental implementations of geomet- 
rically constrained non-linear medium processors, chemical processors and self- 
growing nano-structures. 
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